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Abstract. We report preliminary results of a broadband spectral and temporal study of the black-hole binary 
Cyg X-1 performed with the IBIS telescope. Cyg X-1 was the first pointed celestial target of IBIS during the 
INTEGRAL Performance and Verification Phase, 2002 Nov.-Dec, for a total observing time of ~2 Ms in both 
staring and dithering mode. Here, we report on only the staring, on-axis, observation performed in a stable 
instrument configuration. During the observing period the source was in its characteristic low/hard state, in 
which a few flares and dips have been detected. The IBIS/ISGRI results demonstrate that the INTEGRAL 
observatory offers a unique capability for studying correlations between hardness and/or flux in different bands 
over a wide photon energy range. One of our new results is finding that the hardness-flux correlation changes the 
sign twice over the 20-220 keV; first from positive to negative at ~50 keV, and then back to positive at ~120 
keV. The former change appears to be due to the spectral curvature introduced by variable Compton refiection. 
The latter may be due spectral pivoting. 

Key words, black hole physics - stars: individual: Cyg X-1 - gamma rays: observations - X-rays; binaries - X-rays: 
stars - X-rays: general 



^ 1. Introduction The LH state is characterized by a low flux in soft X- 
K> rays and strong hard X-ray and soft 7-ray flux. Its spectral 
^ Cyg X-1 was the first established black-hole binary (BHB; variability on the time scale of days is then weak, while the 
5^ ; Webster & Murdin 1972; Bohon 1972), and it is so far one ~50 keV flux changes by a factor of a few (e.g., Gierlinski 
of 18 known BHBs and one of three persistent ones. In et al. 1997). The X-ray spectrum in the LH state is weU 
the X-ray band, BHBs usually exhibit five distinct spec- described by the sum of an intrinsic power law with the 
tral /temporal states: low/hard (LH), high/soft (HS), very photon spectral index, F ~ 1.5-2, which is then cut off at 
high, intermediate and quiescent (see a recent review by energies ^ 100 kcV, and a Compton reflection (Magdziarz 
McClintock & RcmiUard 2003). Cyg X-1 itself exhibits & Zdziarski 1995) component. The intrinsic spectrum is 
strong variability on all time scales ranging from millisec- due to thermal Comptonization in a plasma with an elec- 
ond to years, and it is mostly in the LH state. During tron temperature of kT ^ 100 keV and a Thompson op- 
transitions between the LH and HS states, both gradual tical depth of r 1 (Gierhhski et al. 1997), similarly to 
and rapid variations in the hardness ratio have been ob- the case of other BHBs in the LH state (e.g., Zdziarski et 
served. al. 1998). 

Recently, Zdziarski et al. (2002, hereafter Z02) pre- 

T~nr~, TT AD , ■ f sented an analysis of CGi?0/BATSE and RXTE/ ASM 

bend ojjprint requests to: A. Bazzano: angelaOrm. iasf.cnr.it ' ' 

* Based on observations with INTEGRAL, an ESA project observations of Cyg X-1. One of their most interesting re- 

with instruments and science data centre funded by ESA suits m the X/7 range is that most of the LH variability on 

member states (especially the PI countries: Denmark, France, time scales of hundred of days and longer can be explained 

Germany, Italy, Switzerland, Spain), Czech Republic and by varying the slope of the overaU spectrum with a pivot 

Poland, and with the participation of Russia and the USA. energy between 20 and 100 keV. The pivoting can be due 
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Fig. 1. The count rates during the rev. 15-18 vs. the IJD 
(INTEGRAL JuHan Day, starting from 2000 January 1, 
00.00 UT = JD 2451544.5, which corresponds to MJD 
= IJD + 51544). The extent of the data corresponding 
to a given INTEGRAL revolution is marked in the top 
panel. The rates are shown in logarithmic scale with the 
same length per decade for each panel in order to enable 
comparison of the variability in different channels. 



to either variable plasma temperature or variable optical 
depth (Z02). Broad-band X/7 observations are crucial in 
providing constraints on the nature of the variability, and, 
in particular, on the role of the soft photons in cooling the 
X-ray emitting plasma and the role of Compton reflec- 
tion. Here, we present results from the IBIS instrument 
(Ubertini et al. 2003) on board of INTEGRAL (Winkler 
et al. 2003). 



I| 2. Observation and data analysis 

Given the INTEGRAL launch date. Sun aspect con- 
straints led to Cyg X-1 being chosen as the first strong 
7-ray target for the initial calibration of the instruments 
during the Performance Verification (PV) Phase. The ob- 
servations of the Cygnus region in general allowed for ver- 
ification of the imaging performance and the Point Spread 
Function of IBIS as well as its cross-calibration with SPI 
and JEM-X. The PV phase consisted of a staring, on-axis, 
observation lasting ~ 850 ks, and ^ 1130 ks of dithering 
observations to test the SPI performance. Also, 68 ks were 
devoted to calibration of the JEM-X with pointings of 1°- 
6° off axis. In addition, contemporaneous observations of 
Cyg X-1 with RXTE were performed (Pottschmidt et al. 
2003). 

The IBIS instrument (Ubertini et al. 2003) is a coded- 
mask imaging telescope with a large field of view (29° x 
29°) based on two detector layers, ISGRI (Lebrun et al. 
2003) and PICsIT (Labanti et al. 2003) operating in the 
ranges 15 keV-1 MeV and 175 keV-10 MeV, respectively. 
It provides fine imaging (12' FWHM) for source identi- 
fication and good spectral sensitivity in both continuum 
and broad lines over its operative range. During the obser- 
vations, IBIS was still undergoing fine-tuning. Therefore, 
we studied only the staring, on-axis, observations per- 
formed by IBIS in a stable configuration, which corre- 
spond to revolutions 15 to 18, and use the consolidated 
data reprocessed by the INTEGRAL Science Data Center 
(Courvoisier et al. 2003). We use the offline scientific anal- 
ysis (OS A) software v. 1.1. Details and procedures on 
the IBIS software developed by the IBIS Consortium are 
described in Goldwurm et al. (2003). We have analyzed 
ISGRI single interaction events for which arrival time, en- 
ergy deposition and interaction pixels are known. Most of 
the time, PICsIT was in its standard configuration that 
provide images at time interval of the order of several min- 
utes (1750-3600 s). 

We produced the countrate, C, lightcurves during 2002 
Nov. 27-Dec. 8 binned by the pointings (science windows, 
hereafter sews) in the ranges of 20-40, 40-60, 60-120, 120- 
220 and 220-350 keV. The typical duration of a sew is 
~ 2000 s, which allows us to study variability on the time 
scales ranging from that time to days. The highest of the 
above energy channels has not been used due to low sta- 
tistical significance of its single sew. 

The lightcurves in the four energy bands are shown in 
Fig. We see that they are rather similar to each other. 
This is confirmed by calculating the fractional rms vari- 
ability (after subtracting the variability due to the mea- 
surement errors), which is found to be 17 ± 1%, 18 ± 1%, 
16 ± 1%, 14 ± 1% for channels with the increasing energy, 
respectively. 

Fig. [3 shows the corresponding evolution of the hard- 
ness ratio in adjacent channels. As expected by the simi- 
larity of the lightcurves, variability of the hardness ratios 
is much weaker, with the values of the fractional rms of 
3.4 ±0.2%, 2.9 ±0.2%, and 9.3 ±0.5% (with increasing en- 



A. Bazzano et al.: IBIS preliminary results on Cygnus X-1 



3 



ergy) . Interestingly, the spectral variability is strongest at 
the highest observed energies, 60 keV. There is small, 
but statistically significant, increase and decrease of the 
first and second hardness with time. Formally, the best 
exponential fits give the e-folding times of 300 ± 20 days 
and 220 ± 10 days, respectively. On the other hand, there 
is a strong drop in the emission ^ 120 keV from rev. 15 
to 16 followed by a slow decrease during rev. 16, and a 
slow increase during rev. 18. A time scale for significant 
changes of that emission is ^ 3 days. 
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Fig. 2. Hardness ratios for pairs of adjacent channels vs. 
time. The time interval of a given INTEGRAL revolution 
is shown in the middle panel. The ratios are shown in 
logarithmic scale with the constant length per decade in 
order to enable direct comparison of the relative spectral 
variability. 



3. Hardness-flux correlations 

Most of the information on the long-term characteristics 
of Cyg X-1 in hard X-rays/soft 7-rays to date comes from 
the BATSE detector. In the LH state, the 20-300 keV 
data show a positive correlation between the ratio of the 
100-300 keV and 20-100 kcV fiuxes and the 100-300 keV 
fiux. On the other hand, the correlations of the opposite 
sign are seen in the range of the ASM detector, 1.5-12 
keV. Together, the data show the presence of a pivot in 
the - 20-100 keV range (Z02). This is also supported 
by the fractional variability as a function of the energy 
which has a minimum in the 20-100 keV range. As pointed 
out in Z02, the pivoting variability appears only on very 
long time scales, of the order of ^ 100 days. In addition, 
there is a second variability pattern, in which the broad- 
band spectra move up and down with a little change of 
the spectral shape. The presence of both patterns is also 
confirmed by analysis of individual pointed observations 
(Z02, Gierlihski et al. 1997). 

The pivoting variability pattern is explained by Z02 as 
caused by a variable luminosity in the soft seed photons 
irradiating the plasma. If the plasma is not dominated 
by e^ pairs, the variable irradiation results in the plasma 
temperature adjusting itself to satisfy the energy balance. 
As a consequence, softer spectra are expected to corre- 
spond to a lower electron temperature. This appears to be 
confirmed by a BeppoSAX observation in 1996 September 
(Frontera et al. 2001), as pointed out by Z02. On the other 
hand, variations of the local accretion rate in an a fiow of 
a constant geometry may be responsible for the variable 
amplitude with the constant spectral shape. In summary, 
Z02 found strong changes in the hardness on long time 
scale, whereas variability was dominated by changing the 
total luminosity on shorter time scales. 

Following the above suggestions, we consider here cor- 
relations between the spectral hardness and count rates. 
Our results are shown in Fig. El showing the three ra- 
tios (Fig. O each vs. their defining count rates. All the 
hardnesses and count rates are correlated at a very high 
statistical significance (see below). We find a new and un- 
expected result that while the 40-120 keV hardness ratio 
decreases with their defining count rates (top panels in 
Fig. 01, consistent with the pivoting variability of Z02, 
the 20-60 keV hardness increases with the correspond- 
ing count rates. (Note that such a behaviour could not be 
found by Z02, who used only two rather broad BATSE 
channels, 20-100 keV and 100-300 keV.) The change of 
the sign of the hardness-flux correlation in the 40-60 keV 
channel cannot be due to pivoting variability, as hardness 
decreases with flux below a pivot and increases above it, 
just opposite to our result. 

The observed behaviour has to be due to a variable 
spectral curvature, with a hump around ~ 50 keV ap- 
pearing in the spectrum at high fluxes. Such a hump is 
very likely to be due to Compton reflection, well known 
to be present in the spectra of Cyg X-1 (e.g., Gierlihski et 
al. 1997; Gilfanov, Churazov & Revnivtsev 1999). A spec- 
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Fig. 3. Hardness ratios for pairs of adjacent channels as functions of the count rates. The ratios and rates are shown 
in logarithmic scales with the same respective length per decade in each panel in order to enable direct comparison of 
the relative variability. The dashed lines correspond to best fit power-law dependencies. 



trum from reflection of an intrinsic BHB LH-state contin- 
uum is very hard below its peak at several tens of keV 
and very soft above it (e.g., Magdziarz & Zdziarski 1995). 
Thus, its addition hardens the total spectrum below the 
peak and softens it at higher energies, entirely consistent 
with our results. A quantitative study of this effect would 
require data at energies < 20 keV, not available to us at 
this time, as well as spectral fits by each sew, which is 
beyond the scope of this Letter. However, we point out 
that Compton reflection is known to be variable in Cyg 
X-1, and in particular it correlates with the X-ray spectral 
index (Zdziarski, Lubihski & Smith 1999; Gilfanov et al. 
1999). 

Going to higher energies, we also see that the hardness- 
flux correlation continues to be negative up to ~ 120 keV, 
but it becomes again positive at ^ 120 keV. This, on the 
other hand, is consistent with the findings of Z02 (see their 
fig. 6) and may indeed be due to spectral pivoting. 



Fig. O also shows the best fit power-law dependencies, 
corresponding to the count rate ratio cx . The values of 
the exponent are /? = 0.030 ±0.004, 0.073 ±0.004, -0.11 ± 
0.01, -0.10±0.01, -0.30±0.02, and 0.22±0.03, for the six 
consecutive correlations. The small relative uncertainties 
demonstrate the high significance of the correlations. This 
is also confirmed by the Spearman rank correlation test, 
which gives relatitvely high correlation coefficients of 0.29, 
0.43, -0.57, -0.43, -0.42, and 0.24. The corresponding 
probabilities of the absence of a correlation are very low, 
with the highest value of 10~^ for the last correlation and 
<C than that for all others. 

4. Spectral evolution 

We clearly see in Fig.Qla number of short fiarcs and dips 
(note that those are not absorption dips). To study in 
detail evolution of the spectra, we have selected the sews 
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corresponding to the two of the flares and one of the dips. 
For reference, we used stable periods containing 9 and 21 
sews during the rev. 16 and 18, respectively. The log of 
the selected periods is given in Tabled 



Table 1. The log of the observations selected for the study 
of spectral evolution. 



Rev. 


sew 


type 


IJD 


exposure 


16 


6-14 


stable 


1064.6263-1064.9927 


25269 s 


16 


26-27 


peak 


1065.4406-1065.5220 


5428 s 


18 


2-21 


stable 


1070.4397-1071.2539 


59816 s 


18 


26-28 


peak 


1071.4167-1071.5389 


8491 s 


18 


(2)22-24 


dip 


1072.5997-1072.7224 


9557 s 



keV, the optical depth of r ^ 0.9, and the seed photon 
temperature constrained only to ^ 1 kcV. Nevertheless, 
the model parameters are not well constrained due to 
the relatively short exposure, which causes the photon 
statistics to be poor at high energies, ^ 150 keV. To 
overcome this problem, we have tried a simpler model, 
namely, a power law. This yielded the spectral indices 
of r = 2.04 ± 0.02, 2.09 ± 0.02, 2.14 ± 0.02 for sew 25, 
26, 27, respectively. Although xt are always worse than 
those of the comptt fits, the above numbers clearly show 
a softening of the spectrum with the increasing flux. Fig. 
21 confirms this results by showing the ratios of the count 
rate spectra at sew 26 and 27 to that at the sew 25. We 
see that the softening takes place mostly at the highest 
energies, ^ 90 keV. 



At this stage of the IBIS programme a detailed and 
accurate deconvolution matrix is still under finalisation 
and a spectral analysis has to be carried out with caution. 
In this work we have analysed the spectral data using a 
Crab-calibrated matrix with 64 energy channels. This ma- 
trix represents the Crab data well for integration times not 
exceeding 100 ks and energies ^ 500 keV. We have added 
a 5% systematic error to each PHA channel of the spectra 
analyzed below. Spectral fitting was performed with the 
standard XSPEC v. 11.2 tools. 
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Fig. 4. Spectral evolution during the flare in rev. 18 as 
shown by the ratios of the count rate spectra during sew 
26 (red) and 27 (blue) to that from the sew 25 preceding 
the flare (see Table Pi . 



To study the short term spectral evolution of Cyg 
X-1, we have selected three spectra at the highest ob- 
served peak and just before it. This corresponds to the 
day ^1071.5 during rev. 18 and sew 25, 26 and 27 (the last 
one at the maximum of the flare), see Fig. and Tabled 
The three data sets have been fitted to look for spectral 
evolution on the time scale of ~ 3000 s over a wide en- 
ergy range. A good fit {xt ~ 0.9-1.1) has been obtained 
with a thermal Comptonization model (comptt of Hua & 
Titarchuk 1995) at the electron temperature of kT ~ 50 
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Fig. 5. The average physical spectrum for sew 1-33 of 
rev. 18, fitted by a thermal Comptonization model. The 
bottom panel shows the fit residuals. 



Given the above problems with photon statistics, we 
have obtained a spectrum from an uninterrupted period 
of 105 ks during rev. 18 (sews 1-33, IJD 1070.4019- 
1071.7310) during which the source did not show sig- 
nificant spectral variation. A power law model does not 
provide any good fit, and we used the comptt model, as 
above. The obtained parameters are kT = 64 ± 9 keV, 
T = 0.7 ± 0.4 and the seed photon temperature con- 
strained to ^ 1 keV, which yields xt — 1-5 with 51 
d.o.f.. For comparison, the Crab fit with a power law yields 
r = 2.17 ± 0.01 at the 1 kcV normalization of 11.5 ± 0.5 
cm^^ s^-^ at xt — 1-3 for 56 d.o.f. The average physical 
spectrum is shown in Fig. 
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Fig. 6. Count rate ratios between different spectra listed 
in Table 1. Top: The ratio between spectra during the sta- 
ble period in rev. 18 to that in rev. 16. Middle: The ratio 
between the peak in rev. 18 to that in rev. 16. Bottom: 
The peak-to-dip ratio during rev. 18. 



To further study correlations between the flux and the 
spectral hardness, we have measured the count-rate ratios 
between the five spectra listed in Table ^ This method 
allows us to eliminate possible uncertainties due to spec- 
tral deconvolution. In Fig. |H| we plot the ratio between 
two stable periods, two peaks, and between a peak and a 
dip. The first two ratios show no obvious trends, with the 



pairs of spectra just slightly differening in the slope. On 
the other hand, we clearly see a significant softening with 
the increasing flux in the peak-to-dip ratio (the bottom 
panel of Fig. EJ. However, none of the chosen spectra has 
allowed us to measure any dependence of the high-energy 
cutoff (that would measure the plasma temperature) on 
the flux. 



5. Conclusions 

We have shown that the long, almost uninterrupted, obser- 
vation of Cyg X-1 performed during the PV phase with the 
IBIS broad-band imager represents an invaluable database 
to deeply understand the nature of Cyg X-1 at high en- 
ergies. The presented IBIS data cover the energy range 
from 20 keV to a fraction of 1 MeV and a period of 12 
days with interruptions of 7 hours every 3 days due to the 
Earth radiation belt passage. 

During that period, the flux varied by a factor of ~3, 
including periods of stable emission, flares and dips. We 
have first studied correlations between the spectral hard- 
ness and the fiux. Over the studied photon energy range, 
the correlation was found to change the sign twice over 
the 20-220 keV; first from positive to negative at ~50 
keV, and then back to positive at ^120 keV. The former 
change appears to be due to the spectral curvature intro- 
duced by variable Compton reflection. The latter may be 
due spectral pivoting. 

The characteristic time scale for changing spectral 
hardness is the shortest, ^^3 days, at the highest ener- 
gies, ^ 100 keV. Most of the spectral variability is, in 
fact, taking place at those energies. The dependence of 
the high-energy emission on the overall spectral proper- 
ties is relatively complex. In addition to the hardness-flux 
correlations mentioned above, we have found a spectral 
softening with the increasing flux during flares. 
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